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Abstract 

Although recent Higgs data from ATLAS and CMS are compatible with a Standard Model (SM) 
signal at 2a level, both experiments see indications for an excess in the diphoton decay channel, 
which points to new physics beyond the SM. Given such a low Higgs mass run ~ 125 GeV, another 
sign indicating the existence of new physics beyond the SM is the vacuum stability problem, i.e., 
the SM Higgs quartic coupling may run to negative values at a scale below the Planck scale. In this 
paper, we study the vacuum stability and enhanced Higgs diphoton decay rate in the Zee-Babu 
model, which was used to generate tiny Majorana neutrino masses at two-loop level. We find 
that it is rather difficult to find overlapping regions allowed by the vacuum stability and diphoton 
enhancement constraints. As a consequence, it is almost inevitable to introduce new ingredients 
into the model, in order to resolve these two issues simultaneously. 
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I. INTRODUCTION 



In the Standard Model (SM) of particle physics, the Higgs mechanism provides an expla- 
nation to the spontaneous electroweak symmetry breaking, but the Higgs boson itself left 
no trace in all the previous high-energy collider experiments. Recently both ATLAS and 
CMS collaborations announced an observation of a Higgs-like boson at the 5a confidence 



level The discovery points to a new resonance near the 125 GeV, which, if confirmed to 
be the SM Higgs boson, will be a milestone for fundamental particle physics. We assume 
it is the SM Higgs boson in this paper. Such a low Higgs mass then immediately leads to 
problems for vacuum stability, which requires the Higgs self-coupling A remain positive at 
all scales. Given M# ~ 125 GeV, the Higgs self-coupling A would turn negative at a scale 
A ~ 10 9 — 10 11 GeV Q|, indicating the existence of new physics beyond the SM around that 
scale. 

It's quite interesting to notice that both the ATLAS and CMS collaborations report a 
significant enhancement in the diphoton decay channel 

[a(gg — > h) x BR(h — > 77)] 

obs 1 m 1 r\ nn / -i \ 

= [a(gg -+ h) x BR(fc 7 7)]sm = ? °^ ' { ) 

while a(gg — > h) x BR(h — > ZZ*) and a(gg — > h) x BR(/i — > WW*) seem compatible with 

the SM predictions [l|, 0] • In the absence of direct signals of new physics at colliders, the 

enhancement of T(h — > 77) is an important hint of underlying new physics, since the Higgs 

boson couples to photons indirectly via loops where new physics can enter. 

In general, the running of the Higgs self-coupling may be changed by the following types 

of interactions: quartic scalar interaction, Yukawa interaction and gauge interaction. One 

possibility to resolve the SM Higgs vacuum stability problem is to extend the SM with 

additional quartic scalar interactions associated with new scalar degrees of freedom. Typical 

n n 

examples are the Higgs-portal dark matter models, e.g. the darkon |4( and inert |5| dark 
matter models. For a detailed analysis of the implications of these models on Higgs vacuum 



stability, see Refs. |6Ml3j . Another possibility to improve the SM Higgs vacuum stability is 
to introduce new gauge degrees of freedom, in this case one may require new electroweak 
multiplets, the neutral component of which can be cold dark matter candidate, and (or) a 



new U(l) gauge symmetry 13|, in which the SM Higgs carries non-zero charge. To explain 
the enhanced H — > 77 decay width one needs to extend the SM with charged particles fl^ ]. 
which can be charged scalars in the type-II seesaw models , or multi-charged components 



in a septuplet dark matter model 16| that couple to the SM through the quartic scalar 
interaction, or charged vector-like fermion doublets that have Yukawa interactions with the 
SM Higgs, or new W' vector boson which may couple to the SM Higgs through the kinematic 
term or through the W'^W'H^H term. There have been extensive studies focusing on this 



topic 
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22]. 



On the other hand, there are other definitive evidence of new physics beyond the SM, 
including neutrino masses and dark matter. The solar, atmospheric, reactor and acceler- 
ator neutrino experiments have provided convincing evidence that neutrinos are massive 
and lepton flavors are mixed. In order to generate tiny neutrino masses, one may ex- 



23j . Through 



tend the SM with right-handed neutrinos that have large Majorana masses 
Yukawa interactions the right-handed neutrinos couple to the SM lepton doublets, three 
active neutrinos then acquire tiny Majorana masses as given by the Type-I seesaw formula: 
M v = -M D M~ l Ml, where M v is the mass matrix of light active neutrinos, Mq is the 
neutrino Dirac mass matrix linking the left-handed active neutrinos with the right-handed 
neutrinos, M R is the mass matrix of right-handed neutrinos. Actually there are three types 



of tree-level seesaw scenarios 
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24 1 as well as three radiative seesaw scenarios 26l 



2J28 



which may generate Majorana masses for active neutrinos, with the help of dimension-five 
Weinberg operator j^. 

In this paper, we study the Higgs vacuum stability problem and the enhanced H — y 77 
decay rate in a radiative neutrino mass model, i.e. the minimal Zee-Babu model \ 
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which two SU(2) L singlet scalar fields <p + and k ++ are introduced with hypercharge 1 and 
2, respectively. The scalar potential is given as follows 

Vh = (Vh)sm + M 20V + M^K + Ax(0V) 2 + A 2 (/€ t /€) 2 + A 3 (0 t 0)(# t #) 

+X^k)(H^H) + A 5 (0V)(« f «) + + h.c.) , (2) 

where (Vh)sm is the SM Higgs potential. The new Yukawa couplings linking the charged 
scalars to the leptons are 



Cy = f a ^1 a e£ L ^ + + g aP E^ a E Rp K ++ + h.c. 



(3) 



where denote left-handed lepton doublets, E R are right-handed charged leptons. The 
trilinear /1 term in Eq. ([2]) breaks lepton number explicitly, one naturally expects it to be 
reasonably small, since the symmetry is enhanced in the limit fi — > 0. Light neutrino masses 



are generated at two-loop level, which is 



m v = 16nf ac m c g* cd I cd m d f bd , (4) 



where m c are charged lepton masses and iy is a two-loop integral |29| . It turns out that this 
model is not only capable of generating neutrino mass, the interaction terms proportional 
to A3 and A4 also provide a possibility of resolving the SM Higgs vacuum stability problem 
and explaining the enhanced H — > 77 decay rate. 

The paper is organized as follows: In section II, we study numerically the SM Higgs 
vacuum stability and the Zee-Babu scalar potential vacuum stability. Section III is devoted 
to investigating the enhanced H — > 77 decay rate induced by the charged scalar singlets. 
We investigate briefly the searches of the new charged scalars at the LHC in section IV, and 
summarize in section V. 



II. VACUUM STABILITY 

A constraint on the Higgs mass can be obtained by the requirement that spontaneous 
symmetry breaking actually occurs, that is, V(v) is the minimum of the Higgs potential, 
where v ~ 246 GeV is the Higgs vacuum expectation value (VEV) . The one-loop Renor- 
malization Group Equation (RGE) for the quartic coupling, including fermion and gauge 



boson contributions, is 



30] 



levr 2 /^ = +12A 2 - (j<7 2 + 9gVj A + ? (J^gf + \g\gi + g* j + 12\yf - 12yf , (5) 

where the top quark Yukawa coupling is given by y t = s/2M t /v with M t the running mass of 
the top quark. If the initial value of the coupling A is too small, the top quark contribution 
can be dominant and drive it to a negative value A (A) < 0, leading to the scalar potential 
V(A) < V(v). The vacuum is not stable anymore since it has no minimum. Given 
125 GeV, the Higgs self-coupling A would run negative at a scale A ~ 10 9 — 10 11 GeV 

We now analyze the vacuum stability and perturbativity constraints in the Zee-Babu 
model. There are several new parameters, Ai, A2, A3, A4, A5, g and /, which potentially 
affect the RGE running of the Higgs quartic coupling. It turns out that both A 3 and A 4 play 
an important role through their contributions to the RGE of the Higgs quartic coupling. 
These couplings can alleviate both the vacuum stability and perturbativity constraints on 



the SM Higgs mass. In our analysis we employ the two- loop RGEs for the SM couplings [30] 
and one-loop RGEs for the new couplings associated with the Zee-Babu model. The one-loop 
/3-functions of the couplings in the scalar potential can be written as 

16vr 2 /3 A = ^ f^gf + 2 -g\g\ + g\\ - Qy? + 9g^j A + 12A 2 + 2A 2 + 2A 2 - 12yf + 12\y? , 

16vr 2 /3 Al = ^gf - ^ 2 A X + 20A 2 + 2A 2 + A 2 - 2Tr[(fff] + 2Tt[/V]A 1 , 

167T 2 /3 A2 = ^gf - l -fg\\2 + 20A 2 + 2\\ + A 2 - 2Tr[(^) 2 ] + ffi^ , 
27 Q Q 

16vr 2 /3 A3 = - -£ 2 A 3 - -^IA 3 + 6AA 3 + 8AxA 3 + 4A 2 + 2A 4 A 5 + 6y 2 A 3 + Tr[/t/]A 3 , 

16vr 2 /3 A4 = ™g\ - ^A 4 - \g\\, + 6AA 4 + 8A 2 A 4 + 4A 2 + 2A 3 A 5 + 6y 2 A 4 + Tr[^]A 4 , 
432 Q0 

16tt 2 /3 A5 = —g\ - -gi\ + 4A 3 A 4 + 8A1A5 + 8A 2 A 5 + 4A 2 + (Tr[/t/] + Trtfg])\ 5 . (6) 

We assume Tr[/t/], Tr[g^g] <C 1, which is consistent with the requirements of the tiny 
Majorana neutrino masses, such that we can safely neglect them in the numerical analysis. 
and k contribute to the /3-function of gauge coupling gi, which can be written as 

K = (to + x ) 91 ' (7) 

Here the SM gauge couplings are normalized based on SU (5), so the electroweak couplings g, 



g' are related to g x , g 2 via g = g 2 and g' = \/3/%i 3JJ. The couplings aj = gf/4:ir are given 
as a 2 , a 3 ) = (0.01681, 0.03354, 0.1176) at the Z-pole |32j]. The full set of j3— functions 
for gauge couplings, top quark Yukawa coupling; the matching condition between the top 
quark pole mass and running mass, and the one-loop matching condition through which 
Mh is determined by the relation to the running Higgs quartic coupling, can be found in 



Ref. (30[. We set the running mass of the top quark to be M t = 172.5 GeV at the scale of 
Mz in our numerical analysis. 

To ensure vacuum stability, the coefficients of quartic scalar couplings must be non- 
negative in any direction and at any renormalization scale. This leads to the following 
constraints on the scalar couplings: 

• If A 3)4 > 0, the vacuum stability condition is 



A,Ai, 2 >0, A 5 + a/4AiA 2 > . (8) 
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FIG. 1. Running of quartic scalar self-couplings in the Zee-Babu model as a function of the cut-off 
scale, where Ai = 0.1, A2 = 0.1, A3 = 0.2, A4 = 0.4, A5 = —0.2 at the electroweak scale. 



If V A 3j4 < 0, the following additional constraints are needed 



-B 



V / 2AAi" - |A 3 | > 0, ^2\\ 2 - |A 4 | > (for A > or A < & — ^ g [0, 1]), 

B 2 



V"2AAi - |A 3 | > 0, y f 2\X 2 - |A 4 | > 0, C - — < (otherwise) , (9) 



where A = (A 3 - A 4 ) 2 - 2A(A X + A 2 - A 5 ), B = 4AA 4 - 2\j + 2A 3 A 4 - 2AA 5 , C = \j - 2AA X . 
These constraints should be valid at any renormalization scale, and are obtained in analogy 
with the method used in Ref. 33] . In our numerical analysis, we also assume the naive 
perturbativity constraint, namely assume A < 8.2 34| and |Aj| < Air (i = 1 • • • 5) from the 
electroweak scale Aew to the cut-off scale A c . In this analysis, we take the cut-off scale A c 
to be the Planck scale. 

Note that the extra scalars <p, k do not couple to quarks, therefore the running of their 
self-couplings do not receive negative contribution from the top Yukawa coupling. As a 
consequence, they can easily reach the Landau pole if their initial input values at the elec- 
troweak scale are too large. This is indeed what we found in our scan. For small values of 
A 4j 2, we found solutions rendering the vacuum stable at all scales up to the Planck scale. In 
Fig. [1] we plot as an illustrative example the running of the scalar quartic couplings in the 
Zee-Babu model. As the plot shows, after including the contribution of extra scalars 0, k 
to the running, the Higgs quartic coupling A, which turns negative at 10 9 — 10 n GeV in the 
SM, becomes positive up to the Planck scale. From our scan we find that for the vacuum to 
be stable, the couplings A3 and A 4 tend to be of the same sign at the electroweak scale. As 
we will see in the next section when we discuss the enhancement of Higgs diphoton decay 
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rate, this is not the phenomenologically favored case. 



III. ENHANCEMENT OF THE h ->■ 77 DECAY RATE 

As we have mentioned in the introduction, the singly and doubly charged scalars <p 
and k in the Zee-Babu mdoel provide a possible explanation for the enhancement of the 
h — y 77 decay rate. More specifically, <fi and k carry no color charge, implying their van- 
ishing contribution to the gluon fusion production of Higgs, which is the dominate pro- 
duction mechanism at the LHC. However, due to the scalar interaction with Higgs, i.e. 
— \ 3 (<j)^ <p) (H^ H) — \±(k>k)(WH), these two new scalars can augment the Higgs diphoton 
decay rate to some extent. In this section we will show that the large enhancement ratio 
R-y-y recently observed at the LHC can be easily explained in the Zee-Babu model by the 
and k loops. To make our analysis solid and reliable, we will consider the constraints on the 
parameters from perturbativity and collider search results. 



on, k off 



We first turn off the contribution from the scalar k artificially (or this can be understood 
as the limit A4 = at the electroweak scale), and consider the effect of <fi on the Higgs 
diphoton decay rate, to see if is sufficient to produce the large enhancement ratio. The 



partial width for the diphoton channel in the SM is 
T(h -> 77) 



35) 



(10) 



where is the Fermi constant, a the fine-structure constant, Nc = 3 the number of colors, 
Qt = +2/3 the top quark electric charge in units of |e|, r, = 4M 4 2 /M|, i = t,W, and the 
contribution from spin-0, 1/2 and 1 particles are 



A (x) 



A 1/2 (x) 



AAx) 



-x 



2x' 



—x 



X \x 



X \x J \x 



-T + - + 3 - 



x- 



X 



X 



II/- 

X 



(11) 

(12) 
(13) 



respectively. The function f(x) is defined as 

r 

[arcsin a/x] 2 

/(*) = 



f 

4 



In 



i + vT^ 



X' 



1 - VI -x- 1 



for x < 1 



for x > 1 



(14) 



For a 125 GeV Higgs, the contribution from the SM W boson A\ = —8.32 is the dominant 
part and interferes destructively with the contribution from the top quark NcQ 2 Ai/ 2 = 1.84. 

After electroweak symmetry breaking, the scalar acquires mass: M? = Mf + I/2X3V 2 , 
and its coupling to the SM Higgs is g^^x = A31;. With contribution from the new scalar^, 



the enhancement ratio of the Higgs diphoton decay rate is, compared to the SM value 14] 

1 + -^3 



R 



77 



(15) 



2 ■ > M 2 A l (r w ) + N c Q 2 t A 1/2 (r t ) 
It is obvious that the enhancement ratio is determined by two parameters in this case, the 
0-Higgs quartic coupling A3 and the mass M^. We first consider the case A3 < 0, where 
the 0-loop interferes constructively with the SM W boson contribution and tends to enlarge 
the diphoton decay width. 

We should notice that both A3 and are subject to theoretical and experimental con- 
straints. For A3, it is obvious in the left panel of Fig. [2] (A4 = 0) that in some region of the 
parameter space A3 can not be too large, or A becomes non-perturbative at energies much 
lower than the Planck scale. If we consider further the vacuum stability, A3 is more strictly 
constrained^ . For simplicity, we naively assume [ A3 1 < Air. 

In high energy collider experiments, the scalar <fi in the Zee-Babu model behaves like a 
charged Higgs with weak isospin I 3 = whose Yukawa interaction is solely the coupling to 
left-handed leptons. The experiments ALEPH, DELPHI and L3 at LEP have searched for 
charged Higgs 36]. They considered the two or three decay channels H + — > t + z/ t , as, W*A 
{A being light pseudoscalar), giving a limit of 79.3 GeV at 95% confidence level, which is 
independent of the branching ratio BR(ff + — > t + v t ). By analyzing the dilepton events with 
missing transverse momentum, the experiment OPAL obtained a more stringent limit of 
92 GeV at 95% confidence level 3^], assuming the charged Higgs decaying totally to tau 
1 If wc require the potential is bounded from blow, A3 has to satisfy the relation A§ < 2AAi, when the 
scalar k is turned off. Furthermore, even we have Ai > at the electroweak scale, according to the RGE 
running, Ai may become negative at higher energy scales, destabilizing the vacuum. It is a bit non-trivial 
to stabilize the vacuum in the Zee-Babu model. 
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leptons plus missing energy. If the Zee-Babu model is assumed to produce the tiny neutrino 
masses, either for normal hierarchy (NH) or inverted hierarchy (IH) of neutrino spectrum, 
the branching ratios of are determined by the three neutrino mixing angles #^ and the 



Dirac CP violating phase S 



38 



39j, 



NH : ^- 



IH : 



f er 


tan 6*i2 cos #23 


//XT 


COS 013 


fe/i 


tan 612 sin 6*23 


fflT 
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sin 9 23 c _ i5 
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tan #13 


fe/i 


_ COS 023 _« 




tan #13 



+ tan #13 sin#23e tS , 



tan #13 cos#23e" 



(16) 



(17) 



If we utilize the OPAL result to constrain M^, the limit has to be brought a bit down. 
Nevertheless, as the neutrino spectrum is not yet pinned down, we conservatively take the 
orginal OPAL limit of 92 GeV as the lower bound for M$. Charged Higgs has also been 



40]. 



searched for at hadron colliders with higher energies in decays of top quark t — > bH 
However, since the extra scalars <p, k do not couple to quarks, these searches are inapplicable 
for the Zee-Babu model. Therefore, the collider constraint on mass is taken as > 
92 GeV, at 95% confidence level. 

We show in the left panel of Fig. |2]the plots of the enhancement ratios R§ = 1.1, 1.3, 1.5, 
1.7, 2 for h — > 77 rate, as functions of A3 and M^. The blue region on the left is the excluded 
region > 92 GeV. To interfere constructively with the SM contribution, the value of A3 
is negative, as aforementioned. It is clear from the figure that the large ratio _R 77 can be 
easily achieved for a light scalar with mass of hundreds of GeV and a negative coupling 
A3 ~ — 3. 

If A3 > 0, the contribution of the scalar to h — > 77 decay rate interferes destructively 
with that from the SM W boson. To obtain the large h — > 77 width enhancement observed 
at LHC, Mtf, has to be close to but larger than the collider constraint, i.e. 92 GeV, and 
meanwhile the magnitude of A3 is large ~10, i.e. close to its perturbativity bound. Therefore 
the case with positive A3 is not phenomenologically favored. 
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FIG. 2. Left panel: Plots of the enhancement ratios i?^ 7 = 1.1, 1.3, 1.5, 1.7, 1.9, 2 for — > 77) 
as functions of A3 and Ma,, when we turn on only the scalar <f> in the Zee-Babu model, in addition 
to the SM particles. Right panel: contours of = 1.1, 1.3, 1.5, 1.7, 1.9, 2 as functions of A4 and 
M K , when we turn on only the scalar k. The blue regions on the left side in the panels are the 
bounds of collider experiments on Ma, and M K , respectively. See text for details. 

B. 4> off, k on 

In this subsection, we consider the case in which is turned off (or this can be understood 
as A3 = at the electroweak scale), and h — > 77 receives only the contribution from the 
scalar k besides the SM particles. The scalar n carries double electric charges, and the 
corresponding enhancement factor is, analogous to Eq. ([To]) . 

(18) 



K 71 



1 1 v 
1 + -A 4 - 



4A (r K i 



2 l M%A 1 (T W ) + N c Q 2 t A 1/2 (n 

The constraint on A4 from perturbativity is similar to that for A3 discussed in the previous 
subsection, i.e. | A4 1 < Air, regardless of the difference of electric charges carried by 4> and K - 
At colliders, the scalar k acts like "right-handed" doubly charged Higgs with weak isospin 
^3 = ("right-handed" means that the scalar couples to right-handed charged leptons or 



left-handed anti- chargeo 



and hadron colliders 



42 



eptons) . Doubly charged Higgs has been searched for at LEP 41 1 
43[, through the process qq 7*/Z* H ++ H— ^t^^d 



a, b, c, d = e, \i, r), including lepton flavor violating decays, and the strongest bound is 



from CMS at the LHC 



43|. In the inclusive search by CMS, they considered all the six 



decay channels, i.e. BR(H ++ e + e + /^ + /i + /r + r + /e + /i + /e + r + //i + r+) = 100%, as well as 

10 



four benchmark points with different branching ratios, such as 



BR(H ++ -»■ e + e+) = 1/3, 
BR(H ++ -> = 1/3 , 

BR(H ++ -> r + r + ) = 1/3 . (19) 



Based on the pair production analysis only, the exclusion limit on the mass of the doubly 
charged Higgs varies from 165 GeV to 391 GeV, largely depending on the decay mode. In the 
Zee-Babu model, the Yukawa couplings of the doubly charged scalar k to charged leptons, 
g a p, are related to the masses of charged leptons and neutrino oscillation parameters, and 



38 



are constrained by low energy lepton flavor violating processes, e.g. r~ — > /i + yT fi 
The values of the couplings g a p are yet undetermined, we choose the weakest limit on doubly 
charged Higgs from CMS as the mass limit for k: m K > 165 GeV. 

It is worth to mention that the doubly charged Higgs in the CMS analysis is assumed to 



be from the type-II seesaw model 15[ and consequently has weak isospin J3 = 1, which is 



different from the "right-handed" scalar k with I3 = 0. Due to their different coupling to the 



virtual Z boson, the pair production cross section o f rig . 



it-handed doubly charged scalars is 



somewhat smaller than that of the left-handed ones |42l. |44| . When we utilize the mass limit 
originally for the left-handed doubly charged scalar to constrain the mass of right-handed 
ones, we should properly scale down the pair production cross section, thus reduce the mass 
limits. However, to compromise to some extent the dependence of the mass limit on the 
branching ratios of the doubly charged Higgs, we still use the original value of 165 GeV from 
CMS as our limit for m K . 

The plots of the enhancement ratios R*l = 1.1, 1.3, 1.5, 1.7, 2 as functions of A4 and 
m K are shown in the right panel of Fig. [2j The blue region on the left is the constraint 
m K > 165 GeV from CMS. As the electric charge carried by the scalar is doubled, the 
allowed region for the scalar mass and the quartic coupling is somewhat larger, as we expect, 
compared to that for 0. 

If A4 is positive, due to the stringent bound on the scalar mass, as in the case for 0, to 
obtain the large enhancement ratio i? 77 , the scalar mass m K has to be close to but larger than 
its experimental bound, and A4 is very large A4 ~ 10, which is disfavored phenomenologically. 
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C. (j) on ? K on 



Now we turn on both scalar fields (j) and k. In this case, the enhancement ratio is 

2 



1 v 2 A (t^) 1 v 2 AA (t k ) 

1 + ~ A 3^rjy—n T - ; — at 7v> A ? \ + o A 4" 



2 6 MlA 1 (r w ) + N c Q 2 t A 1/2 {r t ) 2 Al(tw) + N C Q 2 A 1 / 2 (r t ) 



(20) 



Collider constraints on and m K and the bound on A3 and A4 from perturbativity remain 
the same. 

We first consider the phenomenologically most favored case in which the two couplings 
A3, 4 are both negative, i.e. —An < A3, 4 < 0, to have contributions interfering constructively 
with the SM W boson. With collider constraints on the scalar masses, > 92 GeV 
and m K > 165 GeV, we now show the plots for the parameter regions allowed to obtain 
the large enhancement ratio i? 77 . We first assume the scalar masses to be close to their 
experimental bounds. More specifically, and m K are, respectively, less than 5 GeV larger 
than their collider lower bounds, denoted as = 92( +5 ) GeV and m K = 165( + 5) GeV. 
We parametrize the two quartic couplings as A3 = (2A4, where a = 5, 3, 1, 1/3, 1/5. The 
enhancement ratios i? 77 as functions of A3 is shown in Fig. [3] as the yellow, red, black, blue 
and green curves. The lighter bands are due to the variations of the scalar masses in the 
regions 92 < < 97 GeV and 165 < m K < 170 GeV. The horizontal dashed black line 
is the enhancement ratio recently observed at the LHC. Obviously, as the two scalars are 
relatively light, IA3 4I are both very small, ~ (9(0.1 — 1). 

When <fi and k are heavier, i.e. with masses of hundreds of GeV, | 4 1 should be somewhat 
large to produce a large _R 77 , as seen in Fig. HJ which shows in the left (right) panel the 
contour of A3 (A4) required to obtain the enhancement ratio i? 77 = 1.71 as a function of 
and m K , with couplings A 3 = 3A 4 (A 4 = 3 A3). When both A 3 , 4 are positive, the situation is a 
bit better than the case with solely one of the scalars, but to produce a large enhancement 
ratio Ryy, (j) , k axe both relatively light and A3, 4 are both large, A3, 4 > 5. This is still not 
the phenomenologically favored case. 

We have to mention that it is almost impossible to find the parameter space where all 
the Higgs vacuum stability, perturbativity and enhanced Higgs diphoton decay rate can be 
achieved simultaneously. To fulfill the Higgs vacuum stability and perturbativity constraints, 
scalar couplings should be small. To achieve the demanded h — > 77 enhancement factor, 
A3 4 should be relatively large. This paradox seems improbable to be resolved without 
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m0=92 (+5) GeV, m /f =165( +5) GeV, A 3 =aA 4 




FIG. 3. Enhancement ratios i? 77 as function of A3 are shown as the yellow, red, black, blue and 
green curves, when we set m$ = 92( +5 ) GeV, m K = 165( + 5) GeV, and parametrize A3 = aX^, with 
a = 5, 3, 1, 1/3, 1/5 in the plot. The lighter bands are due to the variations of the scalar masses. 
The dashed black horizontal line denotes the large enhancement ratio i? 77 recently observed at the 
LHC. 

introducing new ingredients into the model. 



IV. SEARCHES OF <p, k AT THE LHC 



In this section we briefly sketch the search for <fr, k at the LHC. The phenomenology of 
extra scalars in the Zee-Babu model at the LHC has been analyzed in Ref . 39| . Most of the 



discussions in that reference essentially apply here. Recall that the extra scalar fields <fi and 
k do not couple to quarks, the main production of <p and k pairs at the LHC is via quark- 
antiquark annihilation (see Fig. An important feature of this production mechanism is 
that the cross section depends only on one unkonwn parameter, the mass of or k. The 
mass of and k has been constrained to be roughly above 100 GeV by LEP searches for 
charged scalars. In Fig. |6]we plot the LO 4>/k pair production cross section at 8 and 14 TeV 
at the LHC, as a function of their mass. As shown in the figure, the cross section decreases 
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FIG. 6. 4>/k pair production cross section at the LHC. 

quite rapidly with the scalar mass from a few hundred fb for m<^ K = 100 GeV to less 

than 0.1 fb for m^ K = ITeV. The pair production cross section for k is larger than that for 
with the same mass, because k is doubly charged. 

After being pair produced, and k will decay to lighter particles B decays to charged 
lepton and neutrino, whereas k can decay to a pair of </> provided that m K > 2m^, or to a 
pair of charged leptons. Due to the presence of neutrino in the decay product of <p, and also 
due to its smaller production cross section, the discovery potential for is less than that for 
a. A clean and distinct signature for k is provided by its decay to electron/muon pair with 
high invariant mass, where the detection of both electrons and muons is quite efficient. To 
estimate the number of events in these decay channels, we take the optimistic scenario in 



Ref. 39[ where all produced k are assumed to decay to electron/muon pair. In this case, 
we have 207 (3) events for an estimated integrated luminosity 20 fb' 1 at i/i = 8 TeV with 
m K = 200 GeV (m K = 500 GeV), and 8264 (231) events for an integrated luminosity 300 fb' 1 
at v/i = 14 TeV with m K = 200 GeV (m« = 500 GeV). 



V. CONCLUSION 



If the observed ~ 125 GeV boson is, indeed, the Higgs boson, then stability of the SM 
electroweak vacuum up to the Planck scale may require that some BSM degrees of freedom 
become active above the TeV scale. The enhanced Higgs to diphoton decay rate also points 
to new physics beyond the SM. In this paper we studied the Higgs vacuum stability and 
2 We assume these extra scalars are not long-lived, indeed the low energy phenomenology and neutrino data 
imply that k cannot be a long-lived particle 3j|. 
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enhanced Higgs to diphoton decay rate in the Zee-Babu model, which was used to generate 
tiny Majorana neutrino masses at two-loop level. We find that it is rather difficult to find 
overlapping regions allowed by the vacuum stability and diphoton enhancement constraints. 
As a consequence, it is almost inevitable to introduce new ingredients into the model, in 
order to resolve these two issues simultaneously. We also calculated the production cross 
section of new charged scalars at the CERN LHC, which does not change much even if new 
degrees of freedom are introduced into the model. If the model could be made compatible 
with current Higgs diphoton decay rate observation by e.g. introducing new degrees of 
freedom, it would be possible to detect these scalars at the LHC in the near future. 
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